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We study the transition from laminar flow to fully developed turbulence for an inertially-driven 
von Karman flow between two counter-rotating large impellers fltted with curved blades over a 
wide range of Reynolds number (10'^ — 10^). The transition is driven by the destabilisation of the 
azimuthal shear-layer, i.e., Kelvin-Helmholtz instability which exhibits travelling/drifting waves, 
modulated travelling waves and chaos below the emergence of a turbulent spectrum. A local 
quantity — the energy of the velocity fluctuations at a given point — and a global quantity — the 
applied torque — are used to monitor the dynamics. The local quantity defines a critical Reynolds 
number ije^ for the onset of time-dependence in the flow, and an upper threshold/crossover Ret 
for the saturation of the energy cascade. The dimensionless drag coefficient, i.e., the turbulent 
dissipation, reaches a plateau above this finite Ret, as expected for a "Kolmogorov"-like turbulence 
for Re — > oo. Our observations suggest that the transition to turbulence in this closed flow 
is globally supercritical: the energy of the velocity fluctuations can be considered as an order 
parameter characterizing the dynamics from the flrst laminar time-dependence up to the fully 
developed turbulence. Spectral analysis in temporal domain moreover reveals that almost all of 
the fluctuations energy is stored in time-scales one or two orders of magnitude slower than the 
time-scale based on impeller frequency. 
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I. INTRODUCTION 

Hydrodynamic turbulenc e is a key feature for many 
comm on problems (Lesicur, llQQOt iTennekes and Lumlevl . 
I1972D . In a few ideal cases, exact solutions of the Navier- 
Stokes equations are available, based on several assump- 
tions such as auto-similarity, stationari ness, or sym metry 
(for a collection of examples, see ISchlichtingL [T979h . Un- 
fortunately, they are often irrelevant in practice, because 
they are unstable. Two of the simplest examples are 
the centrifugal instability of the Taylor-Couette flow be- 
tween two concentric cylinders, and the Rayleigh-Benard 
convection between two differentially heated plates: once 
the amount of angular momentum or heat is too impor- 
tant to be carried by molecular diffusion, a more efficient 
convective transport arises. Increasing further the con- 
trol parameter in these two examples, secondary bifurca- 
tions occur, leading rapidly to temporal chaos, and/or to 
spatio-temporal chaos, then to turbulence. 

Several approaches have been carried in parallel con- 
cerning developed turbulence, focused on statist i cal pr op- 
erties of flow quantities at small scales ( Frischl . Il995l ) or 
taking into account the persistence of cohe rent structures 
in a more determinis t ic po int of view (|Lesieuij . Il990t 
ITennekes and Ltimlevl . Il972f ). One of the major diffi- 
culty concerning a self-consistent statistical treatment 
of turbulence is indeed that it depends on the flow in 
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which it takes place (for instance wakes, jets or closed 
flows). At finite Re, most turbulent flows could still 
keep in average some geometrical or topological prop- 
erties of the laminar flow (for example the presence of a 
Benard-von Karman street in the wake of a bluff body 
whatever Re), which could still i nfluence its sta t istical 
properties (iLa Porta et all . I2OOII : lOuellette et al. I. I2OO6I: 
IZocchi et all . 11994 ). 

Furthermore, we have recently shown in a von Karman 
flow that a turbulent flow can exhibit multistability, first 
order bifurcations and can eve n keep traces o f its history 
at very high Reynolds number ( Ravelet et al.l . [2004 ). The 
observation of this turbulent bifurcation lead us to study 
the transition from the laminar state to turbulence in 
this inertially driven closed flow. 



A. Overview of the von Karman swirling flow 

1. Instabilities of the von Karman swirling flow between flat 
disks 

The disk flow is an example where exact Navier- 
Stokes solutions are available. The original problem of 
the flow of a viscous fluid over an inflnite rotating flat 
disk has been considered by ,von Karman. (1921) . Ex- 
perimentally, the problem of an infinite disk in an infi- 
nite medium is difficult to address. Addition of a sec - 
ond coaxial disk has be en proposed bvi Batcheloii ( 195ll ) 
and IStewartsonl ( 1953f ). A cylindrical housing to the 
flow can also be added. Instabilities and transitions 
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ha ve been extensively studied in thi s system (for instance 
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The basic principle of this flow is the following: a layer 
of fluid is carried near the disk by viscous friction and is 
thrown outwards by centrifugation. By incompressibility 
of the flow, fluid is pumped toward the centre of the disk. 



Since the review of fzandbergen and Diikstral (119871 ). this 



family of flow is called "von Karman swirling flow". In 
all cases, it deals with the flow between smooth disks, at 
low-Reynolds numbers, enclosed or not into a cylindrical 
container. 



2. The "French washing machine": an inertially-driven, highly 
turbulent von Karman swirling flow. 

Experimentally, the so-called "French washing- 
machine" has been a basis for extensive studies 
of very high- Reynolds numb er turbulence in the 
last decade dBourgoin et al.' 

iDouadv et al., 
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_ _ ■ .Rayelet et al . . _ _ , , 

19961: iTiton and Cadotl . 120031 : IZocchi et all . Il99l . To 
reach a Kolmogorov regime in these studies, a von 
Karman flow is inertially-driven between two disks 
fitted with blades, at a very high- Reynolds number 
(10^ < Re < 10^). Due to the incrtial stirring, very high 
turbulence levels can be reached, with fluctuations up to 
50% of the blades velocity, as we shall see in this article. 

Most of the inertially-driven von Karman setups 
studied in the past dealt with straight blades. Von 
Karman flows with curved-blades-impellers were first de- 
signed b y the VKS-team for dynamo action in liqui d 
sodium (jBourgoin et all . l2002t iMonchaux et al.l . 120071 ). 
With curved blades, the directions of rotation are no 
longer equivalent. One sign of the curvature — i.e., 
with the convex face of the blades forward, direc- 
tion — has been shown to be the most favourable 



to dynamo action (Marie et al.l . 120031 : iMonchaux et al.l . 
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2007"; 'Ravelet et aU 120051 ). The turbulent bifurcation 
plavclct ct al., 2 0041 ) has been obtained with the con- 
cave face of the blades forward, direction (— ). In this 
last work, we discussed about the respective role of the 
turbulent fluctuations and of the changes in the mean- 
flow with increasing the Reynolds number on the multi- 
stability. 
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FIG. 1 (a) Sketch of the experiment. The flow volume be- 
tween the impellers is of height H = 1.8 Rc. (b) Impellers 
used in this article. The disks radius is R = 0.925 Rc and 
they are fitted with 16 curved blades: the two different direc- 
tions of rotation defined here are not equivalent. This model 
of impellers has been used in the VKSl sodium experiment 
(|Bourgoin et all . |2002| ) and is called TM60. 



the measurement techniques in §[IT1 The main data pre- 
sented in this article are gathered by driving our experi- 
ment continuously from laminar to turbulent regimes for 
this peculiar direction of rotation, covering a wide range 
of Reynolds numbers. In ij lllll and §|TV] we characterise 
the basic flow and describe the transition from the lam- 
inar regime to turbulence through quasi-periodicity and 
chaos and explore the construction of the temporal spec- 
trum of velocity fluctuations. The continuity and global 
supercriticality of the transition to turbulence is a main 
result of this article. 

In §|V] we obtain complementary data by comparing 
the two different directions of rotation and the case with 
smooth disk. We show how inertial effects clearly dis- 
criminate both directions of rotation at high-Reynolds 
numbers. 

We then summarize and discuss the main results 
in Wn 



B. Outline of the present article 

Our initial motivation to the present study was thus 
to get an overview of the transition to turbulence and 
to check the range where multistability exists. We first 
describe the experimental setup, the fluid properties and 



II. EXPERIMENTAL SETUP 

A. Dimensions, symmetries and control parameter 

The cylinder radius and height are, respectively, Rc = 
100 mm and He = 500 mm. A sketch of the experiment 
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TABLE I Dynamic viscosity /i (10~^ Pa.s) at various temper- 
atures, density p (kg.m~^) at 20°C and achievable Reynolds 
number range for various mass concentrations C of glycerol 
in water. 



can be found in figure[TJa). We use bladed disks to ensure 
inertial stirring. The impellers consist of 185 mm diame- 
ter stainless-steel disks each fitted with 16 curved blades, 
of curvature radius 50 mm and of height /i = 20 mm (fig- 
ure[T|D). The distance between the inner faces of the disks 
is H = 180 mm, which defines a flow volume of aspect 
ratio H/Rc = 1.8. With the curved blades, the directions 
of rotation are no longer equivalent and we can either ro- 
tate the impellers anticlockwise — with the convex face 
of the blades forward, direction (+) or clockwise (with 
the concave face of the blades forward, direction (— ). 

The impellers are driven by two independent brushless 
1.8 kW motors, with speed servo loop control. The max- 
imal torque they can reach is 11.5 N.m. The motor ro- 
tation frequencies {/i ; can be varied independently 
in the range 1 < / < 15 Hz. Below 1 Hz, the speed 
regulation is not efficient enough, and the dimensional 
quantities are measured with insufficient accuracy. We 
will consider for exact counter-rotating regimes /i = /2 
the imposed speed of the impellers /. 

The experimental setup is thus axisymmetric and sym- 
metric towards rotations of tt around any radial axis 
passing through the centre O (7?.7r-symmetry), and we 
will consider here only T^-jr-symmetric mean solutions, 
though mean flows breaking this symmetry do exist for 
these impellers, at le ast at very high-Reynolds numbers 
([Ravelet et all . 1200 j ). A detailed study of the Reynolds 
number dependence of the "global turbulent bifurca- 
tion" is out of the scope of the present article and will be 
presented elsewhere. Also, since we drive the impellers 
independently, there is always a tiny difference between 
/i and /2 and the 7?.7r-symmetry of the system cannot be 
considered as exact. In the following, we will keep us- 
ing this symmetry — very useful to describe the observed 
patterns- but we will keep in mind that our system is 
only an approximation of a 7?.^ -symmetric system. The 
consequences on the dynamics will be analysed in the 
discussion (ij lVLAj) . 

In the following, all lengths will be expressed in units 
of Rc- We also use cylindrical coordinates {r ; z} and the 
origin is on the axis of the cylinder, and equidistant from 
the two impellers to take advantage of the 7?.7r-symmetry 
(see figure [1^). The time unit is defined with the impeller 



rotation frequency /. The integral Reynolds number Re 
is thus defined as Re — 2TifR^v~^ with v the kinematic 
viscosity of the fluid. 

As in previous works ([Marie and Daviaudl . |2004 
iRavelet et all . l2005l |2004| ). we use water at 20 - 30°C 
as working fluid to get Reynolds numbers in the range 
6.3 X 10^ < Re < 1-2 X 10^. To decrease Re down to 
laminar regimes, i.e., to a few tens, we need a fluid with 
a kinematic viscosity a thousand times greater than that 
of water. We thus use 99%-pure glycerol which kinematic 
viscosity is 0.95 x 10"^ m^.s"! at 20°C (Hodgman, 194;^ 
and should be able to study the range 50 < i?e < 900. To 
cover a wide range of Reynolds numbers and match these 
two extreme ranges, we then use different mixes of glyc- 
erol and water, at temperatures between 15 and 35°C. 
The physical properties of these mixtures are given in 
table m where C is the mass percentage of glycerol in the 
mixture. Solutions samples are controlled in a Couette 
viscometer. 

The temperature of the working fluid is measured with 
a platinum thermoresistance (PtlOO) mounted on the 
cylinder wall ({r = 1 ; z = 0}). To control this tem- 
perature, thermoregulated water circulates in two heat 
exchangers placed behind the impellers. Plexiglas disks 
can be mounted between the impellers and heat exchang- 
ers to reduce the drag of the impellers-back-side flows. 
They are at typically 50 mm from the impellers back 
sides. However, these disks reduce the thermal coupling: 
they are used in turbulent water flows and taken away at 
low Reynolds number. 



B. Experimental tools, dimensionless measured quantities 
and experimental errors 

Several techniques have been used in parallel: flow 
visualisations with light sheets and air bubbles, torque 
measurements and velocity measurements. 

Flow visualisations are made in vertical planes illumi- 
nated by approximately 2 mm thick light sheets. We look 
at two different positions with respect to the flow: either 
the central meridian plane where the visualised compo- 
nents are the radial and axial ones or in a plane almost 
tangent to the cylinder wall where the azimuthal compo- 
nent dominates. Tiny air bubbles (less than 1 mm) are 
used as tracing particles. 

Torques are measured as an image of the current con- 
sumption in the motors given by the servo drives and 
have been calibrated by calorimetry. Brushless motors 
are known to generate electromagnetic noise, due to the 
Pulse- Width-Modulation supply. We use armoured ca- 
bles and three-phase sinusoidal output fllters (Schaffncr 
FN5010-8-99), and the motors are enclosed in Faraday 
cages, which enhances the quality of the measurements. 
The minimal torques we measured are above 0.3 N.m, 
and we estimate the error in the measurements to be 
±0.1 N.m. The torques T will be presented in the di- 
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mensionless form: 



A. Basic state at very low Reynolds number 



Kp = T {pRl{2nff)-^ 

Velocity fields are measured by Laser Doppler Ve- 
locimetry (LDV). We use a single component DAN- 
TEC device, with a He-Ne Flowlite Laser (wave length 
632.8 nm) and a BSA57N20 Enhanced Burst Spectrum 
Analyser. The geometry of the experiment allows us to 
measure in one point either the axial component Vz{t) 
or the azimuthal component Ve{t). Though the time- 
averaged velocity field V is not a solution of the Navier- 
Stokes equations, it is a solenoidal vector field, and it is 
axisymmetric. We thus use the incompressibility condi- 
tion V • V = to compute the remaining radial compo- 
nent Vr- 

The measurements of the time-averaged velocity field 
are performed on a {r x z} = 11x17 grid, using a weight- 
ing of velocities by the particles transi t time, to get ri d 
of velocity biases as explained bv iBuchhave et al.l ( 1979h . 
This acquisition mode does not have a constant acqui- 
sition rate, so we use a different method for the acqui- 
sition of well-sampled signals to perform temporal anal- 
ysis at single points. In this so-called dead-time mode, 
we ensure an average data rate of approximately 5 kHz, 
and the Burst Spectrum Analyser takes one sample ev- 
ery single millisecond such that the final data rate is 
1 kHz. For practical reasons, this method is well-suited 
for points close to the cylindrical wall, so we choose the 
point {r — 0.9 ; z = 0} for the measurements in fig- 
ures [21 m and [5l The signals are re-sampled at 300Hz b y 
a "Sample And Hold" algorithm (jBuchhave et al.l . ll979l ). 

Let us now consider the experimental error on the 
Reynolds number value. The speed servo-loop control 
ensures a precision of 0.5% on /, and an absolute pre- 
cision of ±0.002 on the relative difference of the im- 
pellers speeds (/i — /2)/(/i + /2)- The main error on the 
Reynolds number is thus a systematic error that comes 
from the estimation of the viscosity. As far as the vari- 
ation of the viscosity with temperature is about 4% for 
1°C and the variation with concentration is about 5% 
for 1% of mass concentration, we estimate the absolute 
error on Re to ±10% (the temperature is known within 
1°C). However, the experimental reproducibility of the 
Reynolds number is much higher than ±10%. In the 
range 100 < Re < 500 we are able to impose Re within 
±5. 



III. FROM ORDER TO TURBULENCE: DESCRIPTION 
OF THE REGIMES 

This section describes the evolution of the flow from 
the laminar regime to the fully developed turbulence, i.e., 
for 30 < i?e < 1.2 x 10^. This wide-range study has 
been carried for the negative sense of rotation (— ) of the 
propellers. 



At very low-Reynolds number, the basic laminar flow 
respects the symmetries of the problem. It is stationary, 
axisymmetric and 7?.^-symmetric. This state is stable at 
Re ~ 90, where we present a flow visualisation in fig- 
ure [2fa-b). In figure [2fa), the light sheet passes through 
the axis of the cylinder. The visualised velocities are the 
radial and axial components. The poloidal part of the 
flow consists of two toric recirculation cells, with axial 
pumping directed to the impellers. 

The flow is also made of two counter-rotating cells, 
separated by an azimuthal flat shear-layer, which can be 
seen in figure [2jb) where the light sheet is quasi-tangent 
to the cylinder wall. Both the azimuthal and axial com- 
ponent vanishes in the plane z — Q which is consistent 
with the axisymmetry and the 7?.^-symmetry. This flat 
shear-layer is sketched in figure [2le) . A LDV velocity 
field is presented in section |V] (figure [TUb-d) . 



B. First instability - stationary bifurcation 

The first instability for this flow has been determined 
by visualisation and occurs at Re = 175 ± 5 for both di- 
rections of rotation. The bifurcation is supercritical, non- 
hysteretic, and leads to a stationary regime, with an az- 
imuthal modulation of m = 2 wave number. We present 
a visualisation of this secondary state in figure [2](c) , at 
Re = 185. The axisymmetry is broken: one can see the 
TO = 2 modulation of the shear-layer, also sketched in fig- 
ure [^IJe). One can also note that 7?.7r-symmetry is partly 
broken: the bifurcated flow is T^jr-symmetric with respect 
to two orthogonal radial axis only. This first instabil- 
i ty is very s imilar to the Kelvin-Helmholtz instability. 
dNoreeEaD, iool) made a proper theoretical extension 
of the Kelvin-Helmholtz instability in a cylinder. Their 
model is based on the use of local shear-layer thicknesses 
and Reynolds numbers to take into account the radial 
variations in the cylindrical case. 

We observe this to = 2 shear-layer to rotate very slowly 
in a given direction with a period 7500/"^. This cor- 
responds to a very low frequency, always smaller than 
the maximum measured dissymmetry of the speed servo 
loop control between both independent motors fS III.Bp . 
This is probably the limit of the symmetry of our sys- 
tem, i.e., the pattern is at rest in the slowly rotating 
frame where both frequencies are strictly equal (see dis- 
cussion in ij lVLAp . For convenience, we will describe the 
dynamics in this frame. 

The laminar to = 2 stationary shear-layer pattern 
is observed up to typically Re ~ 300 where time- 
dependence arises. 

To investigate the time-dependent regimes, we now 
also perform precise velocity measurements at a given 
point in the shear-layer. We measure the azimuthal com- 
ponent V0 at {r = 0.9 ; z = 0}, using the dead-time 
acquisition mode (see §111]). 



5 




(b) (c) (d) 

FIG. 2 Visualisation and schematics of the basic laminar flow for impellers rotating in direction (— ). The lightning is made 
with a vertical light sheet. Pictures are integrated over 1/25 s with a video camera, and small air bubbles are used as tracers. 
Picture height is H ~2h = lARc. Laminar axisymmetric flow at Re = 90, meridian view (a). Views in a plane near the cylinder 
wall at Re = 90 (b), Re = 185 (c) and Re — 345 (d). The development of the first m — 2 instabilities — steady undulation 

(c) and rotating vortices (d) — is clearly visible on the shape of the shear-layer. We give sketches of the shear-layer for these 
Reynolds numbers in (e). 



Below, we describe and illustrate the observed dynam- 
ics and the building-up of the chaotic and turbulent spec- 
tra. The next section HVl is complementary: we quantita- 
tively characterize the transitions as much as we can, we 
discuss the mechanisms and we finally propose a global 
supercritical view of the transition to turbulence. 



C. From drifting patterns to chaos 

We present time series of the velocity and power spec- 
tral densities at five Reynolds numbers in figure [3l i?e = 
330 ± 5, i?e = 380 ± 5, i?e = 399 ± 5, i?e = 408 ± 5 and 
Re = 440 ± 5. 



1. Oscillation at impeller frequency 

The point at Re — 330 ± 5 is the first point where a 
clear temporal dynamics is observed: a sharp peak in the 
spectrum (figure [SJj) is present at the impeller rotation 
frequency fa = f — emphasized in the Inset of figure[3fa). 
This oscillation exists for higher Re with the same small 
amplitude: it is too fast to be explicitly visible on the 
long time series of figure [31 but it is responsible for the 
large width of the signal line. 

In comparison, a similar measurement performed at 
Re ~ 260 reveals a flat signal with a very low flat- 



spectrum with just a tiny peak — 1/1000 of the amplitude 
measured at Re = 330 — at = / and we have no data 
in-between to check an evolution. 

On the spectra, we observe the first harmonic, but 
never shows the expected blade frequency 8/ nor a mul- 
tiple. So, it is not clear if it corresponds to the basic fluid 
instability mode or just to a small precessing mode due to 
the misaligning of the impeller axis or even to mechanical 
vibrations transmitted to the fluid through the bearings. 
Since the travelling-wave mode of the next paragraph is 
much stronger and richer in dynamics we will consider 
that the signal at /a = / is a "minor" phenomenon, i.e., 
a perturbation of the steady m — 2 mode. 

In figure[3{a) the mean velocity is not zero, but around 
vg = -1-0.17 during the 600 time units of acquisition, i.e., 
during 600 disks rotations. This value of the velocity 
has no special meaning and depends on the phase be- 
tween the fixed measurement point and the slowly drift- 
ing shear- layer IIII.Bp . The measurement point indeed 
stays on the same side of the shear-layer for this time- 
series but, on much longer time scales, we measure the 
m — 2 shear-layer rotation typical period. 

Further observation on signal and spectrum of fig- 
ure [3^-b reveals some energy at low frequency around 
fa — //30, corresponding to slowly relaxing modulations: 
the slowness of this relaxation is the clear signature of the 
proximity of a critical point. 
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FIG. 3 Temporal signals ve{t) measured by LDV at {r = 0.9 ; ^ = 0} and power spectral densities (PSD), at: (a-b) Re — 330, 
(c-d) Re = 380, (e) Re = 399 , (f) Re = 408 and (g-h) Re = 440. fa is the analysis frequency whereas / is the impellers 
rotation frequency. Inset in (a): zoom over the fast oscillation at frequency /. In (e), a small part of the signal is presented 
with time magnification (x4) and arbitrary shift to highlight the modulation at 6.2/~^. Power spectra are computed by the 
Welch periodogram method twice: with a very long window to catch the slow temporal dynamics and with a shorter window 
to reduce fast scales noise. 



2. Drifting/Travelling Waves 

For 330 < Re < 389 the velocity signal is periodic with 
a low frequency fo- This is illustrated at Re = 380 in 



figure c-d) . The mean velocity is now zero: the shear- 
layer rotates slowly such that the measurement point 
is alternatively in the cell rotating with the upper im- 
peller [vg > 0) and in the cell rotating with the lower 
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impeller {vg < 0). Visualisations confirm that this cor- 
responds to a travelling wave (TW) or a drifting pattern 
and also show that the m = 2 shear-layer is now com- 
posed of two vortices (figure HJi) and thus deserves the 
name "mixing-layer" . Along the equatorial line, one no- 
tice that the parity is broken or the vort ices are tilted 
( Coullet and loossl . 119901 : iKnoblochl . \im(t \. The veloc- 
ity varies between —0.3 ^ vg < 0.3. The drift is still 
slow but one order of magnitude faster than the drift de- 
scribed above for the "steady" m = 2 pattern: one can 
see two periods during 600 time units, i.e., fo = //300 
which is very difficult to resolve by spectral analysis ow- 
ing to the shortness of the signal (c/. caption of figure[3]). 
At Re — 380 (figure [3fc-d), the peak at the rotation fre- 
quency is still present, but starts to spread and becomes 
broadband. The power spectral density at frequencies 
higher than 3/ decreases extremely rapidly to the noise 
level. Let us note that T^Tf-symmetry remains only with 
respect to a pair of orthogonal radial axis which rotates 
with the propagating wave. 



3. Modulated Travelling Waves 

For 389 < Re < 408 the signal reveals quasiperiodic- 
ity, i.e., modulated travelling waves (MTW), shown in 
figure [Ife-f at Re = 399 and 408. The MTW are regular, 
i.e., strictly quasiperiodic below Re = 400 and irregular 
above. The modulation frequency — see magnified (x4) 
piece of signal in figure[3^ — is /m = // (6.2 ± 0.2) what- 
ever Re, even above Re = 400. It is much faster than 
the drift frequency {fn ^ //200) and seems to be re- 
lated to oscillations of the mixing-layer vortex cores on 
the movies. 



4. Chaotic regime 

The upper limit of the regular dynamics is precisely 
and reproducibly Re = 400 and there is no hysteresis. 
From the visualisations, we observe that the m = 2 sym- 
metry is now broken. The mixing-layer vortices, which 
are still globally rotating around the cell in the previous 
direction, also behave more and more erratically with 
increasing Re: their individual dynamics includes excur- 
sions in the opposite direction as well as towards one or 
the other impeller. The velocity signal also looses its 
regularity (see figure [3^-g). 

When this disordered regime is well established, e.g. 
for Re = 440 (figure [3^), it can be described as series of 
almost random and fast jumps from one side to the other 
side of the v = axis. The peaks reached by the veloc- 
ity are now in the range —0.4 ^ vg < 0.4. The spectral 
analysis of the signal at Re — 440 (figure[2h), does not re- 
veal any well-defined frequency peak any more. However, 
a continuum of highly-energetic fiuctuations at low fre- 
quency, below fa = f and down to fa = //lOO, emerges. 
A small bump at the rotation frequency / is still visible, 



and a region of fast fluctuations above the injection fre- 
quency also seems to arise. Although we did not carried 
detailed Poincare analysis or equivalent and cannot char- 
acterize clearly a scenario, we find this transition and this 
regime typical enough to call it "chaos" (see also lIV.Bp . 



D. From chaos to turbulence: building a continuous 
spectrum 

Increasing further the Reynolds number, one obtains 
the situation depicted in figure |4l The time-spectrum is 
now continuous but still evolving. We describe the two 
parts, below and above the impellers frequency /. 



1. Slow/ time scales 

The slow dynamics which has already been described 
at Re — 440 (figure [3^, h) could be thought as depending 
only on the largest spatial scales of the fiow. It is well 
built above Re ~ 10'^ (figure [SId). The mean velocities 
corresponding to each side of the mixing layer are of the 
order of ±0.6 at Re = 1.0 x 10"^ and above (figures 2^, c 
and g). The power spectral density below the injection 
frequency seems to behave with a power-law over 
two decades (see discussion §|VLB]) for all these Reynolds 
numbers (figure [3]). The spectral density saturates below 



2. Fast time scales 

However, the fast time scales, usually interpreted as 
a trace of small spatial scales fluctuations, evolve be- 
tween Re = 1.0 X 10^ and Re = 6.5 x 10^. At the 
former Reynolds number, there are few fast fluctuations 
decaying much faster than f~^^^ (figure HKb)) and the 
intermittent changeovers are easy to identify in the tem- 
poral signal in figure Ufa). If Reynolds number is in- 
creased, the fast (small scales) fiuctuations get bigger 
and bigger and behave as a power law over a growing 
frequency range (figure [3^, f, h). The measured slope 
is of order of —1.55 over 1.5 decade, i.e., 10% lower 
than the classical f~^^^. This value of the exponent 
could be ascribe d to the so-called "bottleneck" effect 
(lFalkovichlll994l) and is compa t ible w ith the values given 
bv lLohse and Muller-Groelind (|l995l ) (-1.56 ± 0.01) for 
a Taylor-micro-scale Reynolds number R\ ~ 100, which 
is an estimation for our flow based on the results of 



IZocchi et all (|l994l) . 



IV. QUANTITATIVE CHARACTERIZATION OF THE 
TRANSITIONS 

The various dynamic states encountered have been de- 
scribed and illustrated in the previous section. Now, 



8 




100 200 300 400 500 600 700 800 



10 

1 

10 



Q 10 
W -1 
Q- 10 

-2 

10 

-3 

10 

-4 

10 



Time (r) 




200 400 600 800 

Time (f"'') 





(e) 


Re=2700 





-3-2-10 1 2 

10 10 10 10 10 10 

f /f 




200 400 600 800 1000 1200 



Time (r) 





3 




10 




2 




10 




1 

10 







Q 


10 


W 


-1 


Q. 


10 




-2 




10 




-3 




10 




-4 




1 




10 




3 




10 




2 




10 




1 

10 







Q 


10 


w 


-1 


Q. 


10 




-2 




10 




-3 




10 




-4 




1 




10 




3 




10 




2 




10 




1 

10 







Q 


10 


W 


-1 


Q. 


10 




-2 




10 




-3 




10 




-4 




1 




10 




3 




10 




2 




10 




1 

10 







Q 


10 




-1 




1 




-2 




10 




-3 




10 




-4 




10 




10 





(b) 


Re=1000 





10 



10 10 

f /f 



10 



10 10 

f /f 



10 



10 10 

f /f 



10 



10 





(d) 


Re=1700 





10 



10 





(f) 


Re=3800 





10 



10 







(h) 


Re=6500 






-3 -2 -1 





1 


10 10 


10 


10 1C 



f /f 

a 



FIG. 4 (a-b): Temporal signal vg{t) measured by LDV at {r = 0.9 ; 2: = 0} and power spectral density at Re = 1.0 x 10''. 
(c-d): Temporal signal and power spectral density at Re = 1.7 x lO'^. (e-f): Power spectral densities at 2.7 x 10^ and 3.8 x 10''. 
(g-h): Temporal signal and power spectral density at 6.5 x 10^. Solid lines in the power spectra plots are power-law eye-guides 
of slope —1 and —5/3. Spectra are computed as explained in the caption of figure O 



we wish to analyse some characteristic measurements — 
the amphtude of the velocity fluctuations and their main 
frequencies — , extract thresholds and critical behaviours 
and then address the question of the nature of the re- 
ported transitions. 



A. From order to turbulence: a global supercriticality 

It is known that fully turbulent von Karman flow can 
generate velocity fluctuations of typically 50% of the 
driving impellers velocity. So, we compute the variance 
Vg j,^^ of the LDV-time-series versus the Reynolds num- 
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FIG. 5 Variance of veif) measured at {r = 0.9 ; 2: = 0} 
vs. Re. Solid line: non-linear fit of the form Vg = a x 
{Re - Recf^"^, fitted between Re = 350 and Re = 2500. 
The regression coefficient is R? — 0.990, and the fit gives 
Rec = 328 ± 8 with 95% confidence interval. The intersection 
between this fit and the asymptotic value v1 ~ 0.27 gives 
Ret = 3.3 X 10^ 
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ber. This quantity is homogeneous to a kinetic energy 
and may be referred to as the azimuthal kinetic energy 
fluctuations in the mixing layer. With this method, we 
consider altogether the broadband frequency response of 
the signal. The results are reported in figure [5] for all the 
measurements performed between 260 < Re < 6500. Ex- 
cept at time-dependence threshold, this quantity behaves 
very smoothly: it can be fitted between Re = 350 and 
Re = 2500 with a law in the square root of the distance 
to a threshold Rec — 330 (figure [51): 

rms {Re - Recf/'^. 

Since we will show below that Rec is precisely the 
threshold for time-dependence, we can make here the 
hypothesis that is a global order parameter for 

the transition to turbulence, i.e., for the transition from 
steady flow to turbulent flow taken as a whole. With this 
point of view the transition is globally supercritical. 



B. Transitions from order to chaos 

We now turn to the very flrst steps of the transition 
to time dependence. We monitor the main frequencies of 
the mixing layer dynamics in the TW- and MTW-regimes 
(see ii lIII.C|) . In these regimes, even if only few periods 
are monitored along single time-series, we carefully es- 
timate the period by measuring the time delay between 
crossings of the v = Q axis. These value are reported in 
flgure [6fa) with circles. In an equivalent way, the peri- 
odicity of the travelling of the mixing-layer vortices on 
the visualisations give complementary data, represented 
by squares on the same flgure. 



FIG. 6 (a) Low-frequency fo of the quasi-periodic regime of 
velocity ve{t) measured at {r = 0.9 ; z = Q\ (circles) and 
drift frequency of the m = 2 shear-layer pattern from fiow 
visualisations (squares with high horizontal error bars due to 
poorer temperature control). The solid line is a linear fit of fo 
between the two thresholds Rerw = 330 and Rechaos ~ 400, 
indicated by vertical dotted lines, (b) Zoom of figure [5] The 
dashed line is a linear fit of the lowest data between Re — 350 
and Re — 450. Close to the threshold, it crosses the dash- 
dotted line which corresponds to the velocity due to the drift 
and estimates the level of imperfection. 



1. Onset of time-dependence 

The drift frequency fo of the travelling waves behaves 
linearly with Re above a threshold Rexw very close to 
330. Both measurement methods agree even if the visual- 
isations deserves large error bars in Re at least due to the 
shortness of our records and to a poorer thermal control. 
So, the flt is made on velocity data only. We observe 
some level of imperfection in the quasi-periodic bifurca- 
tion, due to the pre-existing slow drift below ReTw- we 
always observe the mixing-layer to start rotating in the 
sense of the initial drift. 

We show in flgure [Sfb) a zoom of flgure O i.e., the 
amplitude of the kinetic energy fluctuations. We observe 
that both quadratic amplitude flt and linear frequency fit 
converge to exactly the same threshold Rexw = R^c = 
328. We can conclude that the low- frequency mode at 
fa — Id bifurcates at Re — 330 ± 5 through a zero- 
frequency bifurcation for fo. 

The question is thus how the amplitude precisely be- 
haves at onset. There is obviously a lack of data in the 
narrow range 300 ^ Re < 350 (flgure [DJb)). It is due 
to the high temperature dependence of the viscosity in 
this regime (Reynolds varied quite fast even with thermal 
control) and to some data loss at the time of experimen- 
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tal runs. Despite this lack, we present these observations 
because of the consistency of the different types of data 
— visuahsations, LDV, torques — over the wide Reynolds 
number range. The horizontal line Vg — 0.02 in fig- 
ure [6Ub) corresponds to an amplitude of typically 0.15 
for V0, which is produced just by the initial shear-layer 
drift (see the maximum speed in figure [3^) . This value 
is in good agreement with a linear extrapolation over the 
lower range of figure [6l[b) and thus again with an imper- 
fect bifurcation due to the drift. If we reduce the drift by 
better motor frequencies matching, the onset value of vg 
will depend on the longitude of the probe location and 
the parabola of figure [6ljb) could perhaps be observed on 
the m = 2 shear-layer nodes. 



2. Transition to chaos 

The transition to chaos is very sharply observed for 
Re > Rcchaos = 400. There is no hysteresis. Just above 
the chaotic threshold in the MTW regime (figure [5]:), 
the signal sometimes exhibits a few almost-quasi-periodic 
oscillations still allowing us to measure a characteristic 
frequency. The measured values have been also plotted 
on figure [S^a) and are clearly above the linear fit. This 
could reveal a vanishing time scale, i.e., a precursor for 
the very sharp positive/negative jumps of vg reported in 
the chaotic and turbulent regimes. 

We do not clearly observe any evidence of mode lock- 
ing between the present frequencies which are in the pro- 
gression //200 — » //6.2 f and there is no trace of 
sub-harmonic cascade on any of each. This could be 
linked to a three-f requency scenario d la Ruelle-Takens 
(lMannevillel . ll990D . 
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FIG. 7 Dimensionless torque Kp vs. Re in a log-log scale 
for the negative sense of rotation (— ) of the impellers. The 
main data (o) corresponds to the symmetric (s)-flow regime 
described in this part of the article. For completeness, the 
high-torque br anch (★) for Re > 10* corresponds to the (b)- 
flow regime (iRavelet et al.l . [20041 ). i.e., to the "turbulent bi- 
furcation" fsee ^ IIV.D|) . Since both motors do not deliver the 
same torque in this 7^,r-symmetry broken (b)-flow, the aver- 
age of both values is plotted. Relative error on Re is ±10% ; 
absolute error of ±0.1 N.m on the torque. Rec and Ret are 
the transition values computed from the fits of figure [S] The 
single points, displayed at Re = 5 x lO'', correspond to mea- 
surements in water, where Kp is extracted from a fit of the 
dimensional torque i n a ± & x for 2 x 10^ Re 9 x 10^ 
(|Ravelet et al.1 . |2005| ). 



the minimmn of the Kp{Re) curve (figure [7]). One can 
propose that below this Reynolds number, the power in- 
jected at the impeller rotation frequency mainly excites 
low frequencies belonging to the "chaotic" spectrum, 
whereas above Re ^ 1000 it also drives the high frequen- 
cies through the Kolmogorov- Richardson energy cascade. 



Complementary to the local velocity data, information 
can be collected on spatially integrated energetic data, 

i.e., on torque measurements Kp{Re) (figure E]). The 3. Inertial turbulence 



low- Reynolds viscous part will be described below (§ |V| 
as well as the high-torque bifurcated branch fij llV.Dp . In 
the high-Reynolds number regimes, the torque reaches 
an absolute minimum for Re ~ 1000 and becomes inde- 
pendent of Re above 3300. 



2. From chaos to turbulence 

Is there a way to quantitatively characterize the tran- 
sition or the crossover between chaos and turbulence? It 
seems to be no evidence of any special sign to discrimi- 
nate between the two regimes. An empirical criterion we 
could propose would be the completeness of the [fa/ f)~^ 
low- frequency part of the spectrum, clearly achieved for 
Re = 1000 (figure IDd). This region also corresponds to 



The {Re-Recfl'^ behaviour can be fitted through the 
quasi-periodic and chaotic regimes, up to Re ^ 3000. 
Here, the azimuthal kinetic energy fluctuations level sat- 
urates at v'g ^^j, ~ 0.27, i.e., fluctuations of velocities at 
this point of the mixing-layer are of the order of 50% of 
the impeller tip speed. This saturation is also revealed by 
the Probability Density Functions (PDF) of vq presented 
in figure [H These PDF are computed for 16 Reynolds 
numbers in the range 2.5 x 10*^ ^ < 6.5 x 10^. 
One can notice the bimodal character of the PDF: the 
two bumps, which are symmetric, correspond to the two 
counter-rotating cells. Furthermore, all these PDF col- 
lapse and are therefore almost independent of Re in this 
range. This is also consistent with the spectral data of 
figure mb-d) where the {fa/f)~^ slowest time-scales re- 
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FIG. 8 Probability density function (PDF) of ve for 16 
Reynolds numbers in the range 2.5 x lO'* < Re < 6.5 x 10^. 



gions which contain most of the energy — below / — ap- 
pear similar for Re = 1.0 x 10'^ and above (figure H]). 
The crossover Reynolds number Ret at which the kinetic 
energy of fluctuations saturates in figure [5] is estimated 
by taking the intersection of the horizontal asymptote 
with the fit: Ret — 3.3 x lO'^. This value corresponds 
precisely to the value where the asymptotic plateau is 
reached in the Kp vs. Re diagram (figure [7]). In such 
an inertially-driven turbulent flow, the bulk dissipation 
is much stronger than the dissipation in boundary layers 
and the global dimensionless quantities thus do not de- 
pend on the Reynolds number past a turb ulent threshold 
(|Cadot et all 1 19971: iLathrop et al] . ll992D . 



D. Higher-Reynolds number: multistability and turbulent 
bifurcation 

From all the above reported observations in the nega- 
tive direction of rotation, we conclude that the transition 
to turbulence is completed at Ret and that the azimuthal 
kinetic energy fluctuation can be considered as an order 
parameter for the global transition, from the onset of 
time-dependence Rbc = 330 to the fully turbulent state 
transition/crossover at Ret = 3.3x10^, i.e., over a decade 
in Reynolds number. 

In the inertial regime above Ret, the von Karman flow 
driven by high-curvature bladed impellers rotating in the 
negative direc tion presents another original behaviour: 
iRavelet et al.l (J004) have shown that the turbulent von 
Karman flow can exhibit multistability at high-Reynolds 
number. To study and analyse this phenomenon, it is 
necessary to introduce an additional parameter with re- 
spect to the present paper study: the rotation velocity 
difference A/ = /2 — /i between the two impellers. The 
so-called "Turbulent bifurcation" and multistability are 
observed exclusively for the negative direction of rota- 
tion. So, the A/ = regime presented along this paper 



— called (s) for symmetric in lRavelet et al. can 
be observed only if both motors are started together, i.e., 
if A/ is kept equal to zero at any time. Once some ve- 
locity difference is applied long enough — depending of 
the magnitude of |A/| — , the flow changes abruptly to a 
one cell flow with axial pumping towards one of the im- 
pellers only instead of towards each impe ller. This n ew 
flow — called (b) for bifurcated in .Ravclct et al.l ( 20041 ) — 
strongly breaks the 7?.,r-symmetry, has no middle shear- 
layer and requires much higher torque from the motors: 
typically 3 times the value of (s)-flow, with a finite differ- 
ence between the two motors. The mean reduced torque 
at A/ = is plotted with stars in figure [T] branches 
(s) and (b) co-exist for Re > Re^ = 10^. To recover 
the 7?,7r-symmetric flow, one should stop the motors or at 
least decrease Re below Rem- 
it is worth noting that this multistability is only ob- 
served above Ret, i.e., for flows with a well developed tur- 
bulent inertial Kolmogorov cascade. Furthermore, cycles 
in the parameter plane {Kp2 — Kpi; /2 — /i} h ave been 
made for various Re between 100 and 3 x 10^ (j Raveled . 
|2005| ). At low-Reynolds numbers —Re < 800—, this cy- 
cle is reduced to a continuous, monotonic and reversible 
line in the parameter plane. The flrst apparition of "topo- 
logical" transformations of this simple line into multi- 
ples discontinuous branches of a more complex cycle is 
reported at i?e ~ 5 x 10"^, in the neighbourhood of the 
transitional Reynolds number Ret , and multistability for 
A/ = is flrst observed for Re ~ 10*. The extensive 
study of this turbulent bifurcation with varying Re is 
worth a complete article and will be reported elsewhere. 

From the above preliminary report of our results, 
we emphasize the fact that the turbulent bifurcation 
seems really specific of fully developed turbulent flows. 
Whereas the exa ct counter-rot a ting fl ow (A/ = 0) will 
never bifurcate ([Ravelet et al. I, I2004D . for a smaU A/ 
(0 < |A/|// ^ 1) this turbulent bifurcation around 
Rem — lO'* will correspond to a first order transition on 
the way to infinite Reynolds number dynamics: this flow 
really appears as an ideal prototype of an ideal system 
undergoing a succession of well-deflned transitions on the 
way from order to high-Reynolds-number turbulence. 



E. The regimes: a summary 

The next section concerns some aspects speciflc to the 
inertial stirring. Thereafter is the discussion (ij lVIj) about 
the role of the symmetries and of the spatial scales of the 
flow which can be read almost independently. The fol- 
lowing summary of the observed regimes and transitions 
is given as a support for the discussion. 

• Re < 175 : m = 0, axisymmetric, 7?.7r-symmetric 
steady basic flow (i; IIII.A[) . 

• 175 < Re < 330 : m = 2, discretely 72.7r-symmetric 
steady flow f^ HlLBt . 
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• 330 < i?e < 389 : m = 2, non 

7^7r -symmetric, equatorial-parity-broken travelling 
waves (§ IIII.C.2I § lIV.B.lll , 

• 389 < Re < 400 : modulated travelling waves 

(§inrc3l, 

• 400 < Re < 408 : chaotic modulated travelling 
waves fS IIII.C.3l) . 



• 4 00 < Re < 1000 : chaotic flow (§ IIII.C.4[ 

• 1000 < Re < 3300 : transition to turbulence 

(§ini;Di § iivir2|) , 

• ^ 3300 : inertially-driven fully turbulent flow 

• Re > 10"^ : multivalued inertial turbulence regimes 

(§iiy:d|. 



V. VISCOUS STIRRING VS. INERTIAL STIRRING 

We now focus on the specificities of the inertial stir- 
ring. In the preceding parts, a single rotation sense, the 
negative (— ), was studied. However, very relevant infor- 
mation can be obtained from the comparison of data in 
both senses of impellers rotation, which is equivalent to 
have two sets of impellers with opposite curvature at any 
time in the same experiment. 

The guideline for this analysis is the global energetic 
measurements along the whole Reynolds number range. 
The data for sense (— ) have already been partly dis- 
cussed in the preceding part (figure [Tj), but the full set 
comes here in figure |9l At low Reynolds number the 
two curves are identical, which means that the blades 
have no effect on the viscous stirring. This is analysed 
in § IV. Al However, at high Reynolds number, there is 
a factor 3 between both curves, denoting very different 
inertial regimes, as discussed in § IV.Bl 



the same within 15% for both directions of rotation in the 
laminar regime, wher eas at very high Re , they strongly 
differ (by a factor 2) (jRavelet et al.l . |2005[ ). 

The flow is thus not sensitive to the shape of the im- 
peller blades in the laminar regime. To explain this, we 
make the hypothesis that for these large impellers of ra- 
dius 0.925i?c, fitted with blades of height h = 0.2i?c, 
the flow at low Re is equivalent to the flow between flat 
disks in an effecti ve aspect ratio F = {H — 2h)/Rc = 1.4. 
iNore et all (|2004 !) numerically studied the flow between 
counter-rotating smooth flat disks enclosed in a cylinder 
and report the dependence of the first unstable mode 
wave number on the aspect ratio F = H/Rc. In their 
computations, the critical wave number is m = 1 for 
F — 1.8, whereas for F = 1.4, it is m = 2 as we do 
observe in our experiments. 

We thus compare in figure [TO] our experimental veloc- 
ity fields to a numerical simulation performed by Caroline 
Nore at the same Re and in aspect ratio F = 1.4. The 
three fields are very close. A possible physical explana- 
tion for this effect is the presence of viscous boundary 
layers along the resting cylinder wall. The typical length 
scale of the boundary layer thickness can be estimated 
as S = Re~^/'^. At the Reynolds number when the im- 
pellers blades start to become visible, i.e., at Re ~ 300, 
this boundary layer thickness is of the order of (5 ~ 6 mm, 
while the gap between the impellers and the cylinder wall 
is 7.5 mm. It is also of the order of magnitude of the min- 
imum distance between two blades. For Re < 300, the 
fluid is thus kept between the blades and can not be ex- 
pelled radially: it rotates solidly with the impellers. The 
stirring cannot be considered as inertial and does not 
depend on the blades shape. 

For Re > 300, the dimensionless torque starts to shift 
from a Re~^ law and simultaneously discriminates be- 
tween both sense of rotation: the inertial stirring be- 
comes dominant over the viscous stirring. Simultane- 
ously also, the steady flow becomes unstable with respect 
to time-dependence (§ |TTT] and |TV] ) . 



A. From viscous to inertial stirring 

While Re < 300, the dimensionless t orque Kr, scales as 
Re~^. We are in the laminar regime (ISchhchtind . Il979f ) 
and the viscous terms are dominant in the momentum 
balance. These regimes correspond to ?7i = or m = 2 
steady flows, with an eventual slow drift fi^ lIII.AI fc llll.Bp . 

From the power consumption point of view, both di- 
rections of rotations are equivalent. The two curves — 
circles for direction (— ) and left triangles for direction 
(-I-) — collapse for Re < 300 on a single curve of equa- 
tion Kp = 36.9i?e-^ 

We performed velocity field measurements for the two 
flows at Re ~ 120 - 130 (figure [TOt-f) . The differences 
between the two directions are minor. The order of mag- 
nitude of the mean poloidal and toroidal velocities are 



B. Inertial effects 

At high Reynolds number, we observe in figure [5] dif- 
ferent behaviours for Kp in both sense of rotation. Sense 
(— ) passes a minimum for Re ~ 1000 and then rapidly 
reaches a flat plateau above Ret = 3300 (see § lIV.Cp . 
whereas sense (-f ) asymptotically reaches a regime with 
only a third of the power dissipation of sense (— ). To- 
gether with the curved blade data, figure [9] presents addi- 
tional data for smooth disks. The dimensionless torque 
Kp is approximately 30 times smaller for smooth disks 
than for bladed disks, and does not display any plateau 
at high-Rey nolds number but a Re~^^^ scaling law, as 
described bv lCadot et al. 1 (I1997D . 

It is tempting to compare our curve K „{Re) with the 
classical work of iNikuradsd ([l933, [Ull) consisting in a 
complete and careful experimental data set about the 
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FIG. 9 Compilation of the dimensionless torque Kp vs. Re for various flows. All data stands for 7?.-,r-symmetric von-Karman 
flows except the branch labelled (b— ) (★): see caption of figure[7]for details, (o) : direction of rotation ( — ). (<) : direction of 
rotation (+); the solid line is a non-linear fit of equation Kp = 36.9 x Re~^ between Re = 30 and Re = 250. Some data for flat 
disks of standard machine shop roughness, operated in pure water up to 25Hz (squares) are also displayed with a Re'^^"^ fit. 
Another —1/4 power law is fitted for the positive direction of rotation for 330 < Re < 1500 and is displayed between Re — 10^ 
and Re = 10*. Relative error on Re is ±10% ; absolute error of ±0.1 N.m on the torque. Rec and Ret are the transition values 
computed from the fits of figure (5) 



turbulence in a pipe flow with controlled wall rough- 
ness. The data concern the friction factor — equivalent 
of Kp — measured over a wide range between Re — 500 
and Re = 10^, which is shown to strongly depend of the 
wall roughness above Re ~ 3000. The wall roughness is 
made by controlled sand grains of diameter in the range 
1/507 to 1/15 of the pipe radius, somewhat smaller than 
our blades height h/Rc = 1/5 which can be thought as 
an effective roughness. 

This data set has defied the ory along d e cades and 
still motivates papers. Recen tly, iGoldenfeldl (|2006l ) and 
iGioia and Chakrabortvl ( 20061 ) proposed phenomenologi- 
cal interpretations and empirical reduction of Nikuradse's 
data. In few words, both recent works connect the very 
high-Reynolds inertial behaviour — a plateau at a value 
which scales with the roughness to the power 1 /3 — to the 
Blasius Re~^^^ law for the dissipative region at interme- 
diate Re. IGoldenfeldl ( |2006D . using a method from critical 



point physics, finds a scaling for t he whole domain above 
Re ~ 3000, whereas iGioia andChakrabortv (200^ de- 
scribe the friction factor over the same Reynolds range 
according to Kolmogorov's phenomenological model. 

Compared with pipe flow results and models, our 
Kj,(R e)-cmve (figure El) looks very si milar except for the 
region lGioia and Chakrabortvl (I2006D called the energetic 
regime. Indeed, in our specific case the basic flow itself 
is already dominated by vortices of the size of the vessel. 
The negative direction (circles in figure [5]) shows a min- 
imum followed by a plateau above Ret = 3300 and is in 
agreement with the general inertial behaviour described 
above. However for the positive direction (left triangles 
in figure ini), the Kp curve seems continuously decreasing 
up to Re ~ 10^. Looking closer, one can observe a short 
iT!e~i/'* Blasius regime for Re between 300 and 1500 — 
highlighted by a fit in figure [5] — followed by a very slow 
variation over the next two decades: logarithmic correc- 
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FIG. 10 Comparison betwe en a numer i cal sim ulation (a-b) 
performed with the code of iNore et al.l (|2003l ) in a cylinder 
of aspect ratio F = 1.4 at Re = 120 and two experimental 
velocity fields measured by LDV in direction (+) at Re — 130 
(c-d) and in direction (— ) at Re — 120 (e-f). The fiow quanti- 
ties which we present are in (a-c-e) the azimuthal velocity vo 
and in (b-d-f) the poloidal stream function 'i'. Presenting the 
fields between < r < 1 and < z < 0.7 is sufficient due to 
axisymmetry and 7?,^-symmetry. Blades or smooth disk are 
at 2 = 0.7. 



tions are still visible in the range 10'' < Re ^ 5 x 10^. For 
this direction it is more difficult to d efine a thresh old for 
the plateau observed in pure water ( Mariel . BoOSi ) . Nev- 



ertheless, this threshold should be of order of 10 , i.e., 
much higher than with negative rotation. 

A possible explanation of this strong difference may 
rely in the structure of the flow inside the impellers, i.e., 
in-between the blades. Let us first assume that this flow 
is dominated by what happens along the extrados of the 
blades, on which the pressure is the higher. Then we can 
assume that the blades curvature leads to stable bound- 
ary layers in positive rotation and to Goertler instabil- 
ity in negative rotation. The first case develops Blasius 
boundary layers, whereas the latter develops turbulent 
boundary layers with much more vortices. Therefore, 
when the boundary layer detaches — somewhere along 
the blades or at least at their end — the Blasius boundary 
layer in the positive rotation sheds less turbulent vortices 
than the Goertler's unstable layer does in the negative 
rotation. 

The above description can be sufficient to explain why 
the negative rotation is able to produce a Kolmogorov 
cascade even at quite low- Reynolds numbers near Ret. 
However if, in the positive rotation case, the flow is only 
seeded by vortices produced by the stable boundary layer 
which develops along the smooth blade faces, it is clear 
that a Blasius Re~^^^ can be observed in this transi- 



tion Reynolds range and that a full inertial regime does 
not occur below a very high-Reynolds number owing to 
the very small roughness of the blades faces. This could 
be why both curves in figure [5] look so different: the 
lower one looks qualitatively like a low-roughness bound- 
ary flow and the upper one looks like a high-roughness 
boundary flow. Anyway, this may only account for a part 
of the fiow driving: the resistive torque is much higher 
for any bladed impellers than for flat disks as shown in 
figure [21 

Our observation of the closed von Karma n turbulent 
flow i s thus consistent with the claim by iGoldenfeldl 
(I2006D ■ -hat full understanding of turbulence requires ex- 
plicit accounting for boundary roughness. 



VI. DISCUSSION AND CONCLUSION 
A. Symmetries and first bifurcations 

As for many fiows, the similarity of the fiow behaviour 
at low Reynolds number with intermediate-size non- 
linear system is obvious: breaking a spatial symmetry 
first, then a temporal symmetry and finally transit to 
chaos by a quasi-periodic scenario. 

Comparable study has been carried both experimen- 
tally and numerically in the von Karman flow wi t h fiat 
disk and variable aspect ratio bv lNore et al. I (|2005ll2004 
I2003D . Our results agree well with their results on the 
first instability mode m = 2 if considering the fiuid in 
the blade region as almost solidly driven, which reduces 
the aspect ratio (see IV. Ap . However, all thresholds ap- 
pear at much lower Re for bladed impellers than for fiat 
disks: 175 vs. 300 for the first steady bifurcation and 330 
V.S. above 600 for the first tempor al instability of m = 2 
mode, not observed in lNore et al. I (|2005.1 study. 

Another important difference between both system 
concerns its symmetries. Whereas Nore and collabora- 
tors deal with exact counter rotation by using a single 
motor to drive both disks, our experimental setup uses 
two independent motors and reaches only a approxima- 
tion on a counter-rotating regime. As a consequence, the 
T^-TT-symmetry is stricto sensu broken at any Reynolds 
number and the group of symmetry of our problem is 



SO{2) instead of 0{2). To evaluate the le vel of syinmetr:v 

er (|Chossatl . Il993l 

(/l-/2)/(/l+/2) 



breaking we can use a sma ll parameter (|Chossatl . IIQqJ^ 



iPorter and Knoblochl . [20051) . e.g. e 
which is between 10~^ and 10^^ in our runs. 

Carefully controlling this parameter is an interest- 
ing issue: recently, in almost the same von Karman 
flow in the positive sense of rotation at high Re, 
Ide la Torre and Burguetd ( 2007t ) reported bistabihty and 
a turbulent bifurcation at exactly e = between two 
7^7r-symmetric flows. For non-zero e, the mixing layer 
lies slightly above or below the equator and it randomly 
jumps between these two symmetric positions when e is 
carefully set to zero. 

With our very small experimental e, we verify theo- 
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retica l predictions ( Chossatill993l : |Porter and Knoblochl . 
I2OO5I ) for the 1:2 spatial resonance or fc — 2fc interac- 
tion mechanism with slightly broken reflexion symme- 
try. Instead of mixed mode, pure mode and hetero- 
clinic cycle s — specific of 0(2) and carefully reported by 
iNore et all (|2005l l2004 l2003h — we only observe drift- 
ing instability patterns, i.e., travelling waves and mod- 
ulated travelling waves, characteristic of 50(2). Also, 
the drift frequency is very close to zero at the threshold 
Rcc = 330 (figur e E^), in agreement with the predic- 
tion fo ^ 0{e) (|Chossatl 1 19931: iPorter and Knoblochl . 
I2OO5D . This bifurcation to travelling waves is similar 
to the 1-D drift instability of steady patterns , observed 
in many systems (see, e.g. iFauve et al.L [l99lh . It rehes 
on the breaking of the parity {6 —6) of the pattern 
( Coullet and loosj . ll99Q) : the travelhng- wave pattern is 
a pair of t ilted vortices. The bifurcati on is an imperfect 
pitchfork (|Porter and Knoblochll2005D . 

Finally, the comparison can be extended to the travel- 
ling waves o bserved with flat disks above the mixed and 
pure modes ( Nore et al.ll2005ll2003f ). The observed wave 
frequencies are of the same order of magnitude in both 
case, which let us believe that the same kind of hydro- 
dynamics is involved , i.e., the blades play again a minor 
role at these low Reynolds numbers. However, the fre- 
quency ratio between the basic waves (TW) and their 
modulations (MTW) at onset is much higher 32) in 
our experiment th an in the numerical simulations 5) 
(iNore et al.l . l2003D . This could be due to the high number 
of blades. 

We also wish to consider the symmetry of the von 
Karman flow with respect to the rotation axis. In fact, 
the time-averaged flow is exactly axisymmetric while the 
instantaneous flow is not, because of the presence of 
blades. However, axisymmetry can be considered as an 
effective property at any time at low-Reynolds number 
and at least up to Re — 175, since we have shown that 
the blades have almost no effect on the flow (see § IV. Ap . 
With increasing Re, the blades start playing their role 
and effectively break the axisymmetry of the instanta- 
neous flow. 

Finally, we emphasize that the observations made be- 
low Re ~ 400 closely remind the route to chaos trough 
successive symmetry break for low degree of freedom dy- 
namical systems. Our system can thus be considered as 
a small system — in fact this is coherent with the aspect 
ratio which is of order of 1 — until the Reynolds number 
becomes high enough to excite small dynamical scales in 
the flow. 



B. The three scales of the von Karman flow 

The observations reported in this article — 
visualisations, spectra — evidenced three different 
scales. In particular, time-spectra contain two time- 
frequency domains above and below the injection 
frequency fa = f ■ Let us first make a rough sketch 



of the correspondence between temporal and spatial 
frequency scales of the whole flow: 

• the smallest space-frequencies, at the scale of the 
vessel, describe the basic swirling flow due to the 
impeller and produce the intermediate frequency- 
range, i.e., the peak fa — f in the time- 
spectrum; 

• the intermediate space-frequencies due to the 
shear-layer main instabilities produce the lowest 
time- frequencies; 

• the highest space-frequencies produce, of course, 
the highest temporal frequencies, i.e., the Kol- 
mogorov region. 

The Taylor's hypothesis is based on a linear mapping 
between space- and time-frequencies. It is probably valid 
for the high part of the spectrum, but the mapping might 
be not linear and even not monotonic for the low part. 
We discuss each part of the spectrum in the two following 
paragraphs. 



1. The 1// low-frequency spectrum 

Once chaos is reached at Re — 400, a strong continu- 
ous and monotonic low-frequency spectrum is generated 
(Fig. Oi). In the chaotic regime below Re ~ 1000, the 
spectrum evolves to a neat —1 power law. Then, this 
part of the spectrum does not evolve any more with Re. 

Low-frequency —1 exponents in spectra are common 
and could be due to a variety of physical phenomena: so- 
called "1// noises" have been widely s tudied, e.g., in the 
cond ensed matter field (see for instance iDutta and Hornl . 
I1981D . 

For turbulent von Karman flows driven by two counter- 
rotating impellers, this low time-scale dynamics has been 
alr eady observed over at least a decade in liquid helium 
by IZocchi et al. I (Il994i) as well as for the magnetic in- 
duction spectruni in Hquid metals ( Bourgoin et al.l . l2002l : 
IVolk et al.l . I2OO6I ). However, experiments carried on 
a one-cell flow — without turbulent mixing-layer— did 
not show thi s behaviour (jMariel 120031 : iRaveletl . 120051: 
iRavelet et al. I [2004). We therefore conclude that the 
l//a-spectrum is related to the chaotic wandering of the 
mixing-layer which statistically restores the axisymme- 
try. Once again, the mixing-layer slow dynamics domi- 
nates th e whole dynamics of our sys tem, from momentum 
transfer ( Marie and Daviaudl . l200l to the very high level 
of turbulent fluctuations (Fig. [5] and ^ . 

Furthermore, we can make the hypothesis that the — 1 
slope is due to the distribution of persistence times in 
each side of the bimodal distribution (Fig. [5]) : the low- 
frequency part of the spectrum can be reproduced by a 
random binary signal. Similar ideas for the low-frequency 
spectral construction are proposed for the magnetic in- 
duction in the von Karman sodium (VKS) experiment 
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( Ravelet et al ], l2007l) . In both cases, longer statistics 
would be needed to confirm this idea. 



2. The turbulent fluctuations 

We above emphasize how the flow transits from chaos 
to turbulence between Re ~ 1000 and Ret — 3300. We 
label this region "transition to turbulence" and observe 
the growth of a power-law region in the time-spectra for 
fa > f ■ Does this slope trace back the Kolmogorov cas- 
cade in the space-spectra? 

As the classical Taylor hypothesis cannot apply to our 
full rang e spectrum, we follow the Local Taylor Hypothe- 
sis idea ( Pinton and Labb^ 19941) for the high-fr equency 
part fa > f- Whereas iPinton and Labbd ( 1994 ) did not 
apply their technique — using instantaneous velocity in- 
stead of a constant advection — to the extreme case of 
zero advection, we think it can be applied here owing to 
the shape of the azimuthal velocity PDF (figure[8|) . These 
distributions show first that the instantaneous zero veloc- 
ity is a quite rare event: a local minimum of the curve. 
The modulus of velocity spends typically 75% of the time 
between l/2Vm and 3/2Vm, where ±Vm are the positions 
of the PDF maxima. The sign of the advection has no 
effect on the reconstructed wave number. We can thus 
conclude that frequency and wave number modulus can 
be matched each other at first order by a factor equal to 
the most probable velocity \Vm\ or by the mean of \vg\, 
both very close to each other. This approach is coherent 
with a binary view of the local turbulent signal jump- 
ing randomly between two opposite mean values , just a s 
in turbulent flow reversal model of, e.g., iBenzil HqoI). 
Then, the high-frequency part of the spectrum is equiv- 
alent to the spectrum obtained by averaging the spectra 
of every single time-series between jumps, while the low- 
frequency part is dominated by dynamics of the jumps 
themselves. 

Owing to these arguments, we are convinced that an 
algebraic region dominates the high-frequency part of k- 
spectra above Ret- Observed exponents (—1.55) are of 
the order of the Kolmogorov exponent, less than 10% 
smaller in absolute value. Similar exponents are also en- 
countered at other locations in the vessel. 



C. Conclusion and perspectives 

The von Karman shear-flow with inertial stirring has 
been used for a global study of the transition from order 
to turbulence. The transition scenario is consistent with 
a globally supercritical scenario and this system appears 
as a very powerful table-top prototype for such type of 
study. We have chosen to emphasize a global view over 
a wide range of Reynolds number. This allowed to make 
connections between informations relaying on local (ve- 
locities) or global quantities (torques, flow symmetries), 
as discussed in § IV. Bl and IVIl 



1. Going further 

As a perspective, it would be first interesting to in- 
crease the resolution of the analysis next to the differ- 
ent observed thresholds. It would also be worthwhile to 
perform the same wide-range study for the other sense 
of rotation (-I-) or another couple of impellers. Finally, 
these studies would enable a comparison of the inertial 
effects on the turbulent dynamics at very high Reynolds 
number. 



2. Controlling the mixing layer 

Many results of the present study proceed from veloc- 
ity data collected in the middle of the shear-layer and we 
have shown that this layer and its chaotic/turbulent wan- 
dering can be responsible for the low frequency content 
of the chaotic/turbulent spectrum of the data. 

With the slightly different point of view of controlling 
the disorder level, we have modified the dynamics of the 
shear-layer by addi ng a thin annulus lo cated in the mid- 
plane of the flow ( Ravelet et al.l . l2005t ). This property 
was recently used in the Von Karman Sodium (VKS) ex- 
periment held at Cadarache, France and devoted to the 
experimental study of dynamo action in a turbulent liq- 
uid sodium flow. Dynamo has effectively been observed 
for the first time in this system with a von Karman con- 
figuration using, among other characteris tics, an annulus 
in the mid-plane ( Monchaux et al.l . l2007l ) and is sensitive 
to the presence of this device. Moreover, clear evidence 
has been made that the mixing-layer large-scale patterns 
have a stro ng effect on the magnetic field induction at low 
frequency (jUavelet et all l2007t IVolk et al.l . l2006t) . Fur- 
ther studies of this effect in water experiments are under 
progress. 



3. Statistical properties of the turbulence 

Studies of the von Karman flow currently in progress 
invoke both a wider range of data in space, with the use 
of Stereoscopic 3-components Particle Image Velocimetry 
(SPIV) and a wider range in Reynolds number. 

Whereas the SPIV is slower than LDV and will 
not allow time-spectral analysis, it offers a global 
view of the flow and allows to characterize statistical 
properties of the turbulent velocity. Guided by the 
behaviour of the variance of the local azimuthal velocity 
revealed in the present article (figure O , we expect to 
analyse the evolution of the spatio-temporal statistical 
properties with Re. Such study is very stimulating 
for theoretical advances towar d a statistical mechanics 
of t he turbulence in 2D ( Chavanis and Sommerial . 
Il998l: 



Robert and Sommeria' 



1991) 



(jBouchet an d Sommeria, 

or axisymmet ric flows 

iMonchaux et ah 



200^ 



quasi 



2D 



Jung et al 



(jLeprovost et al 



2006) 



20061 
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